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PREFACE 

This  Memorandum  is  a  product  of  a  study  for  the  Advanced  leccaich 
Projects  Agency  of  Defense  Against  Submarine -Launched  Ballistic 
Missiles.  One  po-tion  of  this  study  is  concerned  with  an  airplane- 
based  boost-intercept  anti-missile  missile  system.  This  Memorandum 
treats  the  problem  of  observation  and  accurate  prediction  of  the 
in-plane  motion  of  an  SLBM  during  the  boost-phase  portion  ol  its 
trajectory,  needed  for  efficient  guidance  of  the  interceptor.  It 
analyzes  the  prediction  accuracy  which  results  when  several  different 
degrees  of  a  priori  knowledge  about  the  boost  trajectory  are  assumed. 
This  Memorandum  is  a  companion  piece  to  an  earlier  Memorandum, 
RM-3606-ARFA,  Early  Estimation  of  SLBM  Heading  for  Boost -Phase .Inter¬ 
ception  (U),  which  considered  the  problem  ol  heading  estimation 
accuracy. 

This  report  should  be  of  interest  to  agencies  and  contractors 
concerned  with  the  interception  of  ballistic  missiles  during  their 
boost  phase. 

The  author  is  a  Professor  in  the  Department  of  Engineering  tuid 
Applied  Science,  Yale  University,  and  is  a  consultant  to  fhe  RAjiD 


Corporation. 
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SUMMARY 

Observations  on  the  trajectory  of  an  SLBM  can  be  used  to  esti¬ 
mate  parameters  of  the  trajectory  and  to  predict  future  positions 
of  the  SLBM  by  using  standard  maximum- likelihood  estimation  methods. 
The  prediction  error  is  a  function  of  the  observation  time  and  the 
length  of  time  for  which  the  prediction  is  made,  and  it  also  depends 
on  the  amount  and  kind  of  a  priori  information  available.  The  pre¬ 
diction  error  is  derived  and  computed  for  four  different  cases  as 
follows: 

1.  The  trajectory  is  a  polynomial  in  time. 

2.  The  trajectory  corresponds  to  constant- thrust  propulsion 
with  unknown  thrust  parameters. 

3.  Constant-thrust  trajectory  with  thrust  parameters  known 
but  time  and  space  origin  unknown. 

4.  Constant- thrust  trajectory  with  only  time  origin  unknown. 

In  general,  the  greater  the  a  priori  knowledge  the  greatei  is  the 
prediction  accuracy. 
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SYMBOLS 

2 

g  =  acceleration  of  gravity  -  ft/sec 

I  *  specific  impulse  of  rocket  fuel  -  seconds 

J  *  covariance  matrix  of  estimated  parameters 

N  =  covariance  matrix  of  noise 

p(s/R)  «*  conditional  probability  density  of  s  given  R 

r  *=  rate  of  fuel  expenditure  -  lbs/sec 

R(t)  =  ground  range  in  ft 

r  =  R(t  )  =  ground  range  at  time  t  *  tc 
c  c 

8R  -  6R(t  )  =  range  error  at  time  t  =  tc 
c  c 

8  =  unknown  time  origin  of  rocket  trajectory 

s  =  true  time  origin  of  rocket  trajectory 
o 

t  *  time  observations  cease  and  prediction  is  made 

i 

t  =  SLBM  booster  cut-off  time 
c 

u  ■  \ 

u  M  «], 

c  4C 

X  *  matrix  of  partial  derivatives  of  range  vs  the  parameters 

W  -  weight  of  SLBM  at  launch  -  lbs 
0 

«  ith  parameter 

a  -  standard  deviation  of  individual  observation  error 
T  -  time  interval  between  observations  -  time  between  radar  pulses 


1 


I .  IOTRODUCTICW 

One  of  the  major  problems  arising  in  the  optimal  guidance  of  an 
interceptor  designed  to  intercept  an  3LBM  during  its  boost-phase  is 
the  estimation  of  the  parameters  of  the  trajectory  of  the  SLBM  and 
the  prediction  of  its  location  at  any  future  time.  1  The  parameters 
can  be  estimated  freer  observations  that  have  been  made  on  the  tra¬ 
jectory,  but  in  general  these  observations  are  not  accurate  because 
of  measurement  noise.  Some  procedure  is  therefore  needed  to  process 
the  data  in  an  optimum  fashion.  The  most  widely  used  method  for  doing 
this  is  the  maximum  likelihood  method,  which  has  been  shown  to  be 
optfum  from  several  points  of  vrev.(2'3)  Therefore,  In  the  following 
analysis  this  method  is  used. 

The  accuracy  of  the  trajectory  prediction  Is  obviously  a  function 
of  the  length  of  time  that  the  SIJ3M  trajectory  has  been  under  obser¬ 
vation  as  well  as  the  amount  of  time  into  the  future  for  which  the 
prediction  is  made.  In  addition  to  these  functional  dependencies, 
the  accuracy  also  depends  very  much  on  the  amount  and  kind  of  a  priori 
information  available  to  the  data  processor.  It  is  this  particular 
dependence  that  is  considered  in  this  Memorandum.  A  number  of 
different  situations  ranging  from  considerable  a  priori  ignorance 
to  almost  complete  a  priori  knowledge  are  analysed. 

That  portion  of  the  trajectory  prediction  problem  which  is  con¬ 
cerned  with  estimating  the  heading  or  azimuth  of  the  trajectory  has 
been  treated  previously. (4)  Since  the  trajectory  is  generally 


nearly  a  straight  line  in  the  latter  portions  of  the  boost-phase,  the 
problem  of  estimating  elevation  is  similar  to  estimating  the  heading. 
In  this  Memorandum  attention  is  focused  on  predicting  the  SLBM 
position  along  its  trajectory  line.  To  simplify  the  discussions,  the 
projection  of  that  notion  on  the  ground  -  the  ground  range  -  vlU  be 
the  quantity  being  predicted. 


II.  PROBLEM 


OTATKMENT  and  assumptions 


The  problem  considered  in  this  Memorandum  is  the  prediction  of 
the  future  ground  range  of  an  JUM  whose  trajectory  is  assumed  to 
lie  in  a  vertical  plane.  It  is  supposed  that  the  range  has  a  mown 
functional  fom  R(t),  and  that  data  are  acquired  and  parameters 
are  estimated  continuously  from  the  time  t  »  0  when  the  missile  is 
first  detected.  The  value  of  H(t)  is  predicted  for  a  future  time 
tc.  The  question  is:  What  is  the  error  in  B(tJ  resulting  from  a 
prediction  made  on  data  acquired  and  processed  up  to  a  time  ^  where 

0  <  ti  <  V 

Although  the  time  t,  sen  have  any  value,  it  will  be  thought  of 
as  the  minimum  cut-off  time  of  the  SUM  booster.  This  might  also 
be  the  intercept  time  under  the  following  conditions: 

a)  The  SLUM  is  on  a  least-favorable  trajectory  .Tom  the  point 
of  view  of  the  interceptor. 

b)  Intercept  is  to  take  place  during  the  boost  phase  of  the 
SLBM  trajectory . 

c)  The  interceptor  guidance  policy  minimises  the  total  fuel 

requirement. 

In  the  examples,  the  following  assumptions  are  made: 

1.  The  only  observed  quantity  considered  by  the  estimator  is 
B(t)  the  ground  range  of  the  missile  from  the  point  of 
first  observation.  This  quantity  is  computed  by  range  and 
angle  measurements  made  by  the  observing  racial , 


2.  The  observations  consist  of  samples  made  at  short,  fixed 

time  intervals.  This  assumption  is  made  primarily  to  simplify 
the  mathematics.  However,  it  also  corresponds  quite  closely 
to  the  operation  of  typical  pulse  radars. 

3.  The  error  in  computed  ground  range  due  to  measurement  error 
at  each  observation  instant  has  a  gaussian  amplitude  distri¬ 
bution  with  zero  mean  and  fixed  variance.  Actually  for 
typical  radars  the  variance  increases  with  distance  from  the 
radar,  but  if  the  difference  between  initial  and  final  dis¬ 
tance  is  relatively  small,  the  change  in  the  variance  is 
also  small. 

4.  The  measurement  errors  from  sample  to  sample  are  uncorrelated. 
This  assumption  is  again  made  for  simplicity.  If  the  measure¬ 
ment  errors  were  assumed  to  be  correlated,  a  form  of  the 

-  correlation  would  have  to  De  sjieclfled.  Also  if  the 'band¬ 
width  of  the  radar  receiver  is  approximately  matched  to  the 
width  of  the  radar  pulses,  there  will  be  very  little 
correlation  between  the  samples. 


III.  RESULTS 


The  following  four  cases  are  considered: 

1.  R(t)  can  be  approximated  by  a  polynomial  in  t  in  which  the 
coefficients  multiplying  powers  of  t  are  unknown,  but  where 
the  tine  origin  is  known. 

2.  R( t)  corresponds  to  a  constant- thrust  trajectory  in  which 
all  parameters  are  unknown.  A  subcase,  in  which  one  of  the 
unknown  parameters  (the  initial  velocity)  is  assumed  to  be 
known,  is  also  considered  under  this  heading  to  simplify 
the  computations  and  to  penult  a  check  with  results  obtalne, 
with  the  digital  computer. 

3.  R(t)  corrcsponus  to  a  constant-thrust  trajectory  in  which 
only  the  initial  value  and  the  tine  origin  are  unknown. 

4.  R(t)  corresponds  to  a  constant-thrust  trajectory  in  which 
only  tnc  time  origin  is  unknown. 

CASE  1.  R(t)  IS  KHOWN  TO  BE  A  CUBIC  POLYNOMIAL  IN  t 

Here  it  is  assumed  that  the  range  function  R(t)  can  be  approxi¬ 
mated  by  a  cubic  polynomial  in  t  during  the  time  interval  0  ■  ■  t  -  t ^ 
Thus 

R( t )  at  ^  a„t“  +  or  t3  (l 

1  C  J 

where  t1 s  are  unknown  parameters.  Typical  trajectories,  such 
as  the  one  assumed  in  Appendix  K  l  Fig-  d)  to  provide  the  common 
basis  for  all  the  numerical  results  of  all  the  analyses  in 
this  Memorandum  con  be  approximated  very  closely  by  a  quadratic 
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polynomial  for  t  <  25  sec,  and  a  cubic  approximation  is  therefore 
good  for  most  of  the  time  interval  of  interest. 

The  computations  for  the  prediction  error  are  straightforward 
and  are  given  in  Appendix  A.  The  ms  prediction  error  is  given  by: 


^Vms 


-  < 


!5<>V 


i  v 


A  graph  of  this  expression  for  tQ  =  50  sec  is  shown  in  Fig.  1.  Note 
that  since  R(t)  is  a  linear  function  of  the  cr's,  the  result  is  inde¬ 
pendent  of  the  o's  (see  discussion  following  Eq.  (?6)  of  Appendix  li). 


CASE  2.  SLBM  MOTION  IS  KNOWN  TO  RESULT  FROM  CONSTANT- THIRST  PRORJLoIOii 
The  constant -thru st  type  of  trajectory  is  more  typical  of  the 
motion  of  SLBM  trajectories  during  the  boost  phase  than  the  poly¬ 
nomial  expression  considered  in  Case  1.  Here  the  acceleration  is 

given  by: 


o 

,  2 

where  g  =  acceleration  of  gravity  in  ft/ see 

I  =  specific  impulse  of  the  rocket  fuel  in  sec 

W  =  initial  weight  of  the  SLBM  in  lb 
o 

r  =  rate  of  fuel  expenditure  in  lb/sec 
This  equation  is  integrated  twi  :c  to  yield 


(3) 
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=  «,  +  o  t  +  1  -  V  log  1  -  \t)  <-  ", 


where  a  =  R(0) .  the  range  at  t  0  in  ft 

at  =  R(0),  the  velocity  at  t  =  O'  in  ft/sec 
=  gl,  ft/sec 
"4  =  i'/w0,  sec'1 

t  =  0  corresponds  to  time  of  first  observation. 

In  this  case  the  ry's  are  all  assumed  to  be  unknown.  The  maximum 
likelihood  estimation  of  the  parameter  and  the  computation  for  the 
predicted  error  are  now  somewhat  more  involved  than  in  Case  1,  since 
R(t)  is  no  longer  linear  in  the  nr's.  Specifically,  the  error  now 
depends  on  the  assumed  true  value  of  04 ,  since  R(t)  is  nonlinear  in 
this  parameter.  (See  discussion  following  Eq.  (po)  of  Appendix  13.) 

The  details  of  the  computations  arc  contained  in  Appendix  3  where 

two  subcases  have  been  considered.  In  the  first  ease  it  was  assumed 

* 

that  the  parameter  ay  -  0.  This  reduces  the  problem  to  on  estimation 
of  three  parameters  and  permits  an  approximate  analytical  solution 
to  be  obtained,  wnicn  holds  for  small  values  of  t  . 

For  large  and  for  the  complete  four- parameter  ’problem  the  com¬ 
putations  were  performed  by  a  digital  computer.  The  composite  results 
are  shown  in  Fig.  <1  for  'y  .01*i3  see’1  and  t  -  ')0  sec--  which  are  the 
parameters  of  the  typical  trajectory  used  in  this  Memorandum  (Fig.  1  of 
Appendix  E) .  The  curve  shows  a  screwhat  larger  prediction  error  than 

"corresponds  to  zero  (or  known)  initial  velocity. 
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the  previous  case,  apparently  because  the  constant-thrust  trajectory 
is  initially  not  very  sensitive  to  changes  in  some  of  the  parameters. 
This  appears  to  result  in  inaccurate  estimates  of  these  parameters. 

CASE  3.  SLBM  MOTION  IS  KNOWN  TO  RESULT  FROM  CONSTANT-THRUST  PRO¬ 
PULSION  AND  ALL  PARAMETERS  EXCEPT  TIME  ORIGIN  AND  STARTING  POINT  ARE 
KNOW  .  . 

The  relatively  large  prediction  error  observed  in  Case  1 
arises  from  the  somevhat  unrealistic  assumption  that  nothing  whatever 
is  known  about  the  parameters  of  the  range  function.  In  practice, 
it  is  quite  possible  that  most  of  the  parameters  are  quite  well 
known,  as  a  result  of  technical  intelligence,  or  otherwise. 

This  applies  especially  to  the  parameters  and  0^  which  can  be 
assumed  to  be  exactly  known  if  the  type  of  rocket  employed  in  the 
SLBM  booster  is  known.  On  the  other  hand,  it  has  been  assumed  in 
the  previous  two  cases  that  the  time  origin  of  the  range  function  is 
known  to  the  observer.  In  general  this  will  not  be  true,  since  the 
missile  is  not  usually  detected  immediately  upon  its  emergence. 

Thus,  in  the  present  case,  it  is  assumed  that  the  range  function  is 
given  by 


where  the  parameters  or  ,  and  ar^  are  assumed  to  be  known  exactly 
and  the  parameters  and  s  are  assumed  to  be  unknown.  In  order  to 
consider  the  possibility  that  some  information  about  (the  range 
from  the  point  at  which  the  SLBM  emerges)  might  be  available  from 


. 


„  independent  radar  observation,  it  is  assumed  that  an  a  priori 
probability  distribution  is  M.  for  V  Tbe  computations  of  this 
case  are  presented  in  Appendix  D,  uhere  it  Is  shove  that  the  a  priori 
probability  distribution  Cor  ^  has  very  little  effect  on  the  pre¬ 
diction  error.  The  results  for  a  typical  set  of  tsirnmeters  are  given 
ln  Fig.  3,  and  they  shou  that  the  prediction  error  is  considerably 
smaller  than  for  Case  2,  'as  would  be  expected. 

CAJEit  m^TON  IS  KM  TO  BE  0F_THg_C0N£^^ 

ONLY~THETDffi  ORIGIN  BEING  UNKNOWN 

This  ease  is  quite  similar  to  the  previous  one  except  that  ^ 
is  now  also  assumed  to  be  known.  The  computations  for  this  case  are 
given  in  Appendix  C,  and  a  graph  of  the  results  is  shown  in  Fig.  k. 
The  ms  prediction  error  is,  of  course,  somewhat  smaller  than  for 


3ase  3. 


A  comparison  of  all  four  cases  is  given  in  Fig.  -uses  3 
and  J.  are  represented  by  only  a  single  curve,  that  for  which  the 
actual  time  origin  is  zero;  i.e.,  it  coincides  with  the  time  origin 

for  the  observations. 


R(t)  corresponds  to  o  constant-thrust 
trajectory  with  oil  parameters  escept 
time  origin  known  _ _ 


t ,  ( sec ) 


Observations  start  ot  t:0 
SLBM  frojectory  starts  ot  t:s0 
tc  -  time  for  which  Rc  is  predicted 
i,  ^  time  observations  cease 
tc  -  s0  =  50  sec 

u4  :  .0143  sec'1  _  _ 


*o 

_ _ L 

S0 :  -5  sec 

s0  ;  -10  sec 


Fig.  4  — Normalized  rms  prediction  error 
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The  computations  of  the  prediction  error  made  In  Section 
III  ignore  an,  a  priori  infection  that  night  he  avallahle,  and 
they  result  in  an  unreasonably  large  prediction  error  for  smal  p 
The  a  priori  information  can  be  taken  into  account  very  roughly  by 
assuming  that  the  prediction  error  is  given  by  the  maximon-llkellhood 
computations  of  Section  III  whenever  this  is  less  than  the  a 
priori  maximum  error,  and  that  the  error  is  othervlse  equal  to  the 
a  priori  maximum.  The  error  in  prediction  then  behaves  roughly  as 
shown  in  Fig.  6,  and  shows  a  pronounced  threshold.  The  time  of  this 
threshold  is  a  function  of  the  a  priori  maximum  error,  and  of  the 
parameters  affecting  the  maximum  likelihood  estimate.  If,  for  example, 
one  assumes  that  the  prediction  error  cannot  exceed  50,000  ft,  and 
lf  0  1.  5000  ft,  and  T  a  .01  sec,  then  the  threshold  occurs  approxi¬ 
mately  at  the  following  times  for  the  four  cases  considered: 


Case 

1 

2a 

2b 

3 

4 


Threshold  Time 
l8.2  sec 
19.8  sec 
22.4  sec 
1.8  sec 
1  sec 


For  values  of  tp  greater  than  the  threshold  value,  the  error  de 
creases  very  rapidly  in  all  cases. 
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V.  CONCLUSIONS 

The  four  examples  of  prediction  of  SUM  trajectories  illustrate 
the  way  in  which  a  priori  information  about  the  trajectory  affects 
the  prediction  error.  If  the  trajectory  is  in  the  form  of  a_poly- 
nomial,  which  means  that  all  time  derivatives  above  a  certain  order 
are  cero,  then  standard  estimation  theory  is  easily  applied  ana  the 
prediction  error  decreases  very  rapidly  with  the  length  of  obser- 
vation  time.  In  practice,  SUM  trajectories  are  often  of  the  constant- 
thrust  type.  This  means  that  their  functional  form,  is  known,  and 
one  would  therefore  expect  that  the  prediction  error  should  be  about 
the  same  as  for  the  case  where  the  trajectory  is  in  the  form,  of  a 
polynomial.  This  is  approximately  true  as  shown  by  the  results  of 
eases  2a  and  2b;  however,  the  error  in  the  constant-thrust  case  is 
uniformly  higher  than  for  the  polynomial.  The  increased  error  1= 
nrobably  caused  by  the  fact  that  the  range  function  is  initially  not 
very  sensitive  to  one  of  the  paremeters.  When  it  is  assumed  that 
the  rocket  parameters  arc  known,  as  in  Cases  3  and  very  small 

prediction  errors  arc  obtained. 

It  is  clear  that  the  results  obtained  in  this  Memorandum 

arc  based  on  the  assumption  that  the  trajectories  depend  on  a 
small  set  of  unknown  but  fixed  parameter.,  and  that  the  functional 
dependence  is  known.  If  actually  a  different  functional  dependence 
from  the  one  assumed  exists,  or  if  important  parameters  arc  omitted 
free,  the  estimation  and  prediction  process,  the  errors  will  obviously 
be  much  ereater  than  those  obtained  here. 


'Phis  observation  applies  particularly  to  the  prosier:  of  predict, 
the  intercept  point  vhen  staging  of  the  3LBM  is  considered,  me 


in  this  ease  depends  heavily,  of  course,  on  the  mount  o 


priori  information  available. 
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Appendix  A 

PREDICTION  OF  SLBM  IK-PLANE  MOTION  WHEN  Hy  t_' 
IS  KNOWN  TO  BE  A  THIRD- OH PER  POLYNOMIAL 


Suppose  that  R(t)  can  be  approximated  during  the  time  interval 


0  <  t  <  t _  by  a  polynomial  expansion  of  the  form: 


R(  t)  -■  a,  t  +  Qf^.t"  • 


(6) 


There  is  no  difficulty  in  principle  in  consider  inf  hidicr -or  icr  ! CD  - 
nomials,  except  that  the  mathematics  become  more  unwieldy.  The  quali¬ 
tative  results  should  not  be  affected  by  the  oidez  Oi  -.  e  ip]* ox....ation 
in  any  ease.  Note  that  Eq.  (6)  implies  Uvit  the  SLBM  is  launched  it 
the  time  t  0,  and  that  the  location  of  the  launch  point  and  the  time 
of  launch  are  known. 

We  assume  that  the  parameters  a,,  and  are  to  be  estimated 
from  data  taken  or.  R(t)  only.  The  acquisition  of  data  and  the  esti¬ 
mation  of  the  parameters  is  assumed  to  take  place  continuously  iron 
the  time  the  target  is  launched.  Then  the  question  is:  What  is  Jie 
error  in  R(  t  )  resulting  from  a  prediction  made  on  data  acquired  an! 
processed  up  to  a  time  t  where  0  <  ^  <  t  1 


For  a  third-order  fit,  we  have  from  Eq.  (0  that  the  cnor  m 


R( tj  ii 


p  3 

6R(t  )  £w,t  to  t  *■  let  t 


( 7 ) 


where  6a.  is  the  error  in  estimating  a,.  Hie  uncertainly  in  hvt(-) 
is  best  described  by  the  variance  giver,  by: 


_0 


1 

“  ** 

8R(t  ) 

^  2  ^  k  h 

=  tc  (8^)  +  2tJ  bcY1  6%,  +  tc 

2 

(8rv2)  +  26~x  5o'j 

+  2t5  dry  5^  +  V  (&y,)J  (8) 

c  2  3  c  j 


vhere  the  bar  indicates  an  ensemble  average. 

It  is  supposed  that  R(t)  is  sampled  at  short  intervals,  ?,  and 
that  the  measurement  error  has  a  normal  distribution  with  zero  mean 
and  a  covariance  matrix  N.  For  convenience  we  let  t  =  ir  where  n 
is  an  integer;  then  N  is  an  n  x  n  matrix.  If  maximum  likelihood 
estimation  is  employed,  then  the  covariance  matrix  of  the  bo's  is 
given  by ^  ^ 


where  the  usual  vector  notation  is  employed,  and  where,  from  Kq.  (i) 


oR( t)  - 


where  1  =  1,  2,  ... 
For  simplicity 


-  1 

»—  €♦ 

_ 1 

t 

0R(t) 

= 

n 

tc 

=  it 

.  . 

t3 

( s 

\  C Hr 

U 

the  transpose 

It 


It 


OR 

Onr 


(10) 


vc  assume  that  the  measurement  errors  on  R  are 


uncorrelated  from  sample  to -sample,  and  have  the  same  variance  cT. 
(Actually,  the  variance  increases  with  distance,  but  for  a  short 
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time  interval  the  change  of  the  variance  is  snail.)  Then  the 
covariance  matrix  N  reduces  simply  to  where  IQ  is  the  n-diaen 

sional  unit  matrix,  and  therefore  the  general  element  JrE  of  the 
matrix  J  defined  in  Eq.  (9)  becomes  simply 


(11) 


The  partial  derivatives  have  the  form  given  in  Eq.  (10). 


For  example 


J 


12 


n 

£  (iT)(iT)2 

i=l 


4 

n 


4 

T 


40  T 


4 


where  we  have  used  the  fact  that  ^  -  nr  and  where  the  error  In  the 
approximation  ia  of  order  l/n.  For  t  .  .01  tec,  thla  error  beeaaon 
negligible  for  l;  >  1  sec.  The  other  elements  of  J  can  be  obtained 
In  a  similar  manner,  giving: 


<L 

0  T 


(13) 
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The  covariance  matrix  of  the  parameter  errors  is  the  inverse  of  this 


-1  2 
J  =  15  a  t 


20 

+  3 


60 

T¥ 


42 

.5 


60 

3 


192 

.5 


140 

.6 


42 


140 


105 

+7 


(14) 


Substituting  the  values  from  Eq.  (14)  into  Eq.  (0)  gives  the  result: 


J 


When  the  square  root  of  Eq.  (ly)  is  plotted  on  log-log  paper 
the  points  for  t^/t  <  l/2  fall  on  a  straight  line  having  a  slope 
of  -3.69.  For  t  -  50  sec,  the  nns  error  is  therefore  given  very 
closely  by 


5R(t  ) 
c 

a  t 


rms 


y.  09  x  1.0f 


.  - 3 ■ 89 
h 


(lu) 


t,  <  1/2  t  . 

1  '  c 


Tills  approximation  holds  for  t 
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Another  approximate  form  for  Eq.  (15)  Is 


By  comparing  coefficients  of  t^  it  can  be  seen  that  this  expres¬ 
sion  differs  from  Eq.  (15)  by  less  than  1  per  cent  for  any  value  of 

t  /t  .  It  has  therefore  been  used  to  plot  the  curves  shown  in  Figs, 
c'  1 

1  and  5- 


Appendix  B 


PREDICTION  OF  SLBM  IN-PLANE  MOTION  DURING  BOOgf 
UHTTj  PROPULSION  IS" KNOWN  TO  BE  CONSTANT  THRUST 


Suppose  that  the  ground  range  R(t)  of  an  SLBM  is  given  by  a 
known  nonlinear  function  of  a  set  of  parameters  a  =  ...  ^1, 


i.e., 


R(t)  =  t{a,  t) 


(18) 


If  the  function  f(o, t)  is  smoothly  varying  in  v,  then  the  error 

6R  =  6R(t  )  of  the  predicted  value  of  R  at  the  time  of  cut-off, 
c  c 

is  related  to  the  errors  in  estimating  the  a' s  by  the  approximate 
relation: 


dfKt  )  8f(2,t  ) 

6Rc  „  6^  *  ♦ 


Df(or,tc) 


oo. 


+  0(6 rr") 


(19) 


The  higher-order  terms  can  be  neglected  if  either  the  5o^  or  the 
higher  partial  derivatives  are  small  enough.  Whether  this  is 
permissible  in  any  particular  case  obviously  depends  on  f(v,t)  and 
the  observation  error,  and  the  approximation  error  should  always  be 
checked.  For  the  moment  we  assume  that  Eq.  (19)  io  valid.  Then  the 
variance  of  the  prediction  error  is 


where  ~  and  T  are  respectively  the  variances  and  covariances 

5o'^  1  J 

of  the  error  in  estimating  the  ft's. 

It  is  supposed  that  the  <r's  are  estimated  on  the  basis  of 
measured  values  of  R(t),  and  that  the  measurement  error  is  nonnaliy 
distributed  with  zero  mean.  It  is  assumed  that  R(t)  is  sampled  at 
short  intervals  t,  and  that  at  the  time  of  estimation  n  samples  of 
R(t)  have  been  obtained;  thus  ^  -  nr.  Then  the  Joint  distribution 
of  all  the  n  measurement  errors  is  normal,  with  an  n-dlmensional 
variance  matrix  N. 

If  the  n's  are  estimated  by  use  of  maximum  likelihood,  then  the 
estimated  value  of  2  is  given  by  the  solution  of  the  maximum  Ukeli- 

(3,8) 

hood  equations : 


M_1 -  f  (")1  =  o 


(21) 


is  smoothly  varying  it  can  be  expanded  in  a  Taylor  series  about  a 


nominal  value  of  2 >  Hq’  as  Tollows: 


f(«)  =  £(lJ  +  x<2  '  2q)  + 


and  Eq,  (21)  can  be  replaced  by  a  linearized  version: 


X  if1  [R  -  HzJ  -  X(2  “  =  0 


which  has  the  explicit  solution: 


ry  -  ry 

-  -0 


%  -  rt  =  x  n"1  X  X  N-1  R  -  r(y  ) 


—  “  —o 


Note  that  the  matrix  X  of  partial  derivative  is  generally  not  square 
and  it  therefore  docs  not  have  an  inverse.  Hence  the  X'n  cannot  be 


cancelled  in  Eq.  (23).  However,  the  matrix  product  XII  ‘  xj  io  a 

k  x  k  symmetric  matrix  which  is  generally  non-singular  and  therefore 
¥ 

has  an  inverse. 

Suppose  now  that  ^  is  actually  the  correct  value  of  v-  Then 

y  -  a  is  the  error  in  the  estimation  of  v.  Also  f(nr  )  is  then  the 
-o 

correct  value  of  R,  and  R  -  £izQ)  *  thc  >*oise.  The  variance  of 
the  estimation  error  is  then 


_  r  -i-l  _  r  1-1 

(a  -  0 'J2  =  I  N"1  X  XN-1  nn  N_1  X  Ml-1  X 


•If  this  matrix  is  singular,  this  is  un  indication  that  thc 
measurements  are  not  sensitive  to  one  or  more  of  the  porametexs.  In 
this  case  the  ratio  of  estimation  error  to  measurement  error  quite 
properly  becomes  infinite. 


where  the  symmetry  of  N  and  of  om*'s‘;,^on 

transpose  symbols  at  the  right.  However,  m  -  N,  the  noise  covariance 

matrix.  Hence 


<2  -  ^ 


)f  N_1  X 


1-1 


-  J 


-1 


(24) 


It  is  necessary  to  check  the  validity  of  the  linear  approximation; 
this  can  be  done  by  using  the  rms  estimation  error  computed  from 
Eq.  (24)  as  an  indication  of  the  actual  error,  and  then  to  use 
this  to  investigate  the  magnitude  of  the  higher-order  terns  of 
Eqs.  (19)  and  (22).  It  is  then  generally  possible  to  place  bounds 
on  the  magnitude  of  the  noise  variance  terms  for  which  the  approxima¬ 
tion  is  Justified.  This  will  be  done  in  the  specific  application 

given  below. 


For  simplicity  we  shall  again  assume  tnat  the  measurement  errors 


in  R(t)  are  uncorrelated  from  sample  to  sample  and  have  the  same 


variance  a2.  In  this  case  the  noise  covariance  matrix  N  reduces  to 


A  where  I  is  the  n-dimensional  unit  matrix.  Then  the  covariance 
-n  —  n 


matrix  J  of  the  estimation  error  become! 


J  -  4; 


X  X 


so  that  the  typical  element  of  J  becomes : 


1  ^  df(o>iT)  dfCn-jii) 

J-  =  -2  I  “ 


rs 


Orr 


OCT 


(25) 
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If  t  is  small,  then  the  summation  can  be  replaced  by  an  integral  so 
that 


h  dfKt)  dfKt) 

1  |  ~  — r -  dt 


on 


ch 


(20) 


a  t 


1  ,  . .  1 

fact  that  Ai—  di  -d(Ti)  7 


dt 


where  use  has  been  made  of  the 
and  that  nr  =  t^. 

It  should  be  noted  that  if  f(2,t)  is  actually-  a  linear  function 
of  some  of  the  o's,  then  neither  the  covariance  matrix  of  the  esti¬ 
mation  error  nor  the  variance  of  the  prediction  error  will  be  a 
function  of  these  <v's.  Thus,  suppose  that  f(2,t)  is  a  linear  function 

of  7^.  Then 


f (2 ,t )  ( t )  o1  '  g(2,t) 


where  g(^,t)  is  not  a  function  of  Then 


of  (or ,  t ) 


(.nr . 


^(t) 


which  is  not  a  function  of  oy  Thus  the  matrix  X  defined  in  Kq.  (•  1, 
will  not  be  a  function  of  «1  and  neither  will  the  error  covariance 
matrix  J_1  defined  in  Kq.  (2).)-  Again,  since  the  partial  derivatives 
c)f(o,tc)/^1  are  not  functions  of  ^  if  f(o,t)  is  a  linear  function 
of  or  ,  the  prediction  error  given  by  Kq.  (20)  is  also  independent  of 

these  q^. 


t  ■■ 


•c. 


1 


ih'Il 


{ :«_■»  v- 


?ftArof 


EXAMPLE  OF  A  CONSTANT -THRUST  TARGET 


Consider  now  on  SLBM  target  on  a  straight-line  constant -thrust 


trajectory.  The  range  is  given  by: 


w  -  rt  /  v  -  rt  \ 

H(t)  =  Ko  .  Voi  i  git  UU—  log  j  *  1 


Rq  =  range  in  feet  at  tine  t  -  0 

Vc  =  velocity  in  ft/sec  at  time  t  =  0 

g  =  acceleration  of  gravity  in  ft/sec6- 

I  =  specific  impulse  of  the  rocket  fuel  in  sec 

w  =  weight  of  the  target  vehicle  in  lb  at  t  =  0 
o 

r  =  rate  of  fuel  expenditure  in  lb/sec 

If  we  consider  the  parameters  a ,  =  R  ,  a  ,  -  V  ,  a  gl,  and 

1  o  o  3  ' 

'4  =  r/w  ,  then 


f(n,t)  =  0f1  +  Ogt  +  —  (1  -  Qr^t)  log  (1  -  orJ(t)  +  arj.t  {21 


c)f(o,Jt) 


Sf(a,t) 


r  -I 

— t - =  —  (1  -  a.  t)  log  (l  -  a.t)  +  a,  t 


32 


T  -  T  -  x 

d23  “  d32  "  2  3 

0  t  uJ 


J2k  =  J42  =  ~2  5 

0  T  U 


.h 

bit. 

J34  Jj+3  ~  1* 

a  t  u 


*3  O 

l6uJ  -  3u  -  ou  1/,  ,  2  .  3\  ,  , .  \ 
- ^ - -  7  (1  -  3u  ♦  2uJ)  log(l-u) 


log(l-u)  -  ^  (4u^  -  3u2  -  6u) 


|(l-u)2  log2(l-u) 


|(l  +  6u  -  9u"  +  2u^)  log(l-u)  +  j^(6u  +  21uu  -  lbu^) 


where  u  =  ar^ 

The  variances  and  covariances  of  the  estimated  parameters  needed 
in  Eq.  (20)  to  determine  6Rc  are,  according  to  Eq.  (2U) ,  the  elements 
of  the  inverse  matrix  j"1.  The  inversion  of  the  1*  x  4  j  matrix  whose 
elements  are  given  in  Eq.  (30)  is  best  done  by  means  of  a  digital 
computer.  However  an  approximate  idea  of  the  solution  can  be  gained 
by  reducing  the  problem  somewhat  and  approximating  the  log  functions 
in  Eq.  (30)  by  a  power  series.  In  this  connection  it  should  be  noted 
that  if  uc  -  o^tc,  1  -  uc  l/A,  where  A  is  the  macs  ratio  of  the  rock¬ 
et.  Hence  0  <  u  <  1.  For  small  t  the  expansion  for  log(l-u)  there¬ 
fore  can  be  expected  to  converge  quite  well. 

The  dimensionality  of  the  problem  can  be  reduced  by  assuming  that 
the  initial  velocity  of  the  target  vehicle  is  known  to  be  zero.*  This 


♦The  computations  given  below  indicate  that  the  major  component 
of  the  prediction  error  is  contributed  by  errors  in  the  last  two  para¬ 
meters  a j  and  0^.  Hence  removal  of  either  or  av,  does  not  have  much 

effect  on  the  results. 


removes  the  second  tern  from  Eq.  (27)  and  eliminate,  the  elements  J^, 
Jj2,  J23,  J24  from  J.  The  approximate  expressions  for  the  remaining 

elements  of  J,  obtained  by  expandlne  loe(l-n)  into  a  pover  Berles 

about  u  =  0,  become: 


JU  "  2 

a  t 


.3 


3  ^ 

J33  2  \  20  3b  b3 

at' 


l  u2 


J44 


•K  ( i 

~T~  l  20 

a  t  \ 


J13  =  J31 


2 

a  t 


+  IB  +  u2  7^j+  - 


2  3 

u  ^  u  u 

l  +  27  +  Zb 


14  =  °4l 


2 

a  t 


J3'>  =  J43  '  a% 


1  U  u 

l  +  12  +  20  + 


2  o„3 

u  2u 

20  +  24  +  63 


(H) 


The  inverBion  of  the  3  x  3  J  matrix  1b  now  fairly  straightforward. 

The  detailed  computations  show  that  the  matrix  is  very  nearly  singular, 
i.e.,  the  determinant  ’J  consists  of  terns  which  approximately  cancel 
each  other,  leaving  only  a  small  residue.  The  computations  for  the 
inverse  must  therefore  be  fairly  exact.  This  fact  must  be  kept  in  mind 
when  the  larger  4X4  matrix  is  computed  by  means  of  the  digital 
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2  f  16128 

,2  or  L  + - 


[6R  r  = 

C 


(l-uc)log(l-uc)  +  uc 


16128 


log(l-uc)  +  uc 


12 


TT 


336  [(l-uc)log(l-uc)  +  u,  33b[log(l-uc)  ♦  uj 


(33) 


32256  (i-uc)log(l-uc)  +  uj  [log(l-uJ  +  uc 


where  u  =  a,  t  .  Note  that  a?  does  not  appear  in  this  result. 

C  4  c  3 

Equation  (33)  has  been  computed  for  the  following  numerical  example: 

t  =  50  sec 
c 

=  .0143  sec’1 


The  computation  shows  that  the  second,  third,  and  sixth  ..eras  of  Eq. 

7 - T5 

(33)  are  very  much  larger  than  the  other  three  terms.  Thus  5Rc 

J.  ^ 


decreases  with  t|;  in  fact, 


6R 


can  be  approximated  by 


_ _  2 

[BR]1*  =  16128  —At 

Vl 


(2  -  uc)  log(l  -  uj  +  2uc 


(34) 


The  curves  shown  in  Figs.  2  and  5  v^re  plotted  from  computer 
results.  The  curve  for  three  unknown  parameters  coincides  quite 
closely  with  the  approximate  expression,  Eq.  (34),  for  values  of 
u  <  .1;  for  large  values  of  u  the  ras  error  predicted  by  Eq.  (34)  1b 
considerably  greater  than  the  value  obtained  by  the  computer  for  both 
the  three-  and  four-parameter  cases. 
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A  check  on  the  validity  of  linear  approximations  used  in  this 


development  was  undertaken  by  computing  the  ratio  r  of  second-order 


terms  to  linear  terms  in  Eq.  (19)  as  applied  to  the  constant-thrust 


propulsion  example.  It  is  clear  from  Eq.  (29)  that  the  only  non-zero 


second  portials  are: 


3^f(a,t) 


*3  °“4 


=  -  4  iogd 
ak 


or^t)  +  or4t 


d‘‘f(a,t)  2(1“0fl+t)  l°g(l-or4t)  +  2a  -  cr^  t‘ 


Hence  the  ratio  r  is  given  by: 


do-j  t-b'jj  3  4 


f  “  a2  (&o'.)t 


'•/A 

a  *r5°'i 


where  f  1b  shorthand  for  f(of,t).  The  ratio  is  evaluated  by  substi¬ 


tuting  the  rms  vuIucb  of  the  5o^  in  this  expression,  and  if  r  is  less 


than  setae  arbitrary  figure  (say,  .1)  the  linear  approximation  may  be 


taken  to  be  reasonably  accurate. 


This  computation  has  been  performed  for  the  following  values  of 


t  »  50  sec 
c 


or i  =  .0143  sec 


c*3  =  1970  ft/Bec 


(35) 


the  parameters: 


3? 


corresponding  to  the  typical  trajectory  U6ed  throughout  this  Memo¬ 
randum  and  given  in  Appendix  E.  It  should  be  noted  that  the  contri¬ 
bution  of  the  second-order  terns  is  no  longer  independent  of  and 
a  value  must  therefore  be  specified.  A  curve  of  a  / 7  vs  t^  for  r  «  .1 
is  given  in  Fig.  ?.  From  this  curve  it  can  be  seen  that  if,  for 
instance,  a  / 7  =  100  ft  sec^  then  the  second-order  terms  in  Eq.  (19) 
are  less  than  10  per  cent  of  the  first-order  terms  for  >  18  sec. 

The  curve  rises  rather  rapidly  with  t  so  that  for  the  same  value  of 
o  /t  the  approximation  error  is  less  than  1  per  cent  for  >  28  sec. 
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Appendix  C 


PREDICTION  CF  ST.RM  IK-PLANE  POSITION  DURING  BOOST  WHEW  PROPULSION 
IS  KNOWN  TO  BE  CONSTANT-THRUST  AND  WHB<  ONLY 
THE  TIME  ORIGIN  IS  UNKNOWN 


It  is  supposed  that  the  ground,  range  of  the  target  missile  is 


given  by 


R  =  R(t  -  s) 


where  s  is  unknown,  but  where  the  functional  form  of  R  is  known 


exactly.  Discrete  observations  are  made  on  R  starting  at  some  time 


t  =  0  and  ending  at  time  t  .  The  time  of  the  first  observation  tQ 


occuirs  some  seconds  after  the  SLBM  is  launched.  At  time  t^  it  is 


desired  to  estimate  R(t  )  -  R  ;  in  general  t  >  t. .  The  observations 

c  c  c  X 


on  R  are  spaced  t  seconds  apart  and  are  corrupted  by  additive  noise, 


which  is  assumed  to  be  white,  gaussian,  and  stationary,  with  zero  mean 


and  variance  a  . 


by  P.  M.  Woodward. 


bility  of  s,  given  the  observed  values  of  R,  is  given  by; 


p( b/R)  -  k  p(s)  e6^ 


where  k  is  a  normalizing  constant  choBen  such  that 


p(s/R)  ds  =  1 


The  problem  of  estimating  s  is  essentially  the  problem  considered 
u  (9;10)  yoodvard  shows  that  the  a  posteriori  proba- 


p( a)  is  the  a  priori  probability  density  of  s  which  is  assumed  to  be 
uniform  in  the  sequel.  Under  the  assumption  of  white,  gaussian, 


the  noise,  and  t  =  it  is  the  i—  observation  instant.  The  upper 
limit  of  summation  n  is  such  that  t^  =  nT  -  t^,  the  time  at  which  the 
observations  cease  and  the  prediction  is  made. 

For  s  near  s  ,  g(s)  can  be  expanded  in  a  Taylor  series. 

2 

g(s)  =  g(s  )  +  (s  -  s  )  -^|  +  J(s  -  s  )2  + 

&SJ  “'o'  v  o'  os|s=s  i  o  dsHs=s 


From  Eq.  (4l)  it  is  clear  that 


6(8.)  =  -h  y  LR2(ti-so)  +  2r(veo}  n(ti} 


2  a 


i=l 


(43) 


n 

Sil.  ■  *'(■„) --4  i-so> 


i=i 


(44) 


Ml 

ds  s  o 

3  1  =  1 


+  n(tt)  R*(VB0) 


(>»5) 
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Woodward^  has  shown  that  if  — ^  R  (t^-sQ)  ?  >  ^  ^ 


20  W 


-  r  -2 

^  ,  R/(ti-SQ)  1  >  >  1,  then  the  a  posteriori  probability  p(s/R) 


can  be  obtained  approximately  by  ignoring  the  terms  of  Eqs.  (43),  (44) 


and  (45)  containing  n(t).  In  this  case,  the  g'(sQ)  term  is  neglected. 


Also,  s(sQ)  is  not  a  function  of  s  and  can  be  combined  with  the 


normalizing  constant  of  Eq.  (39)-  The  expression  for  the  approximate 


a  posteriori  probability  becanes: 


p(s/R)  sik'e 


--T  ^  R,(ti"So) j  (c*Do}i 


where  k'  is  a  new  normalizing  constant  containing  e  and  also 


the  a  priori  probability  p(s)  which  was  assumed  to  be  uniform  (i-e., 


constant).  Equation  (46)  shows  the  a  posteriori  probability  to  be 


approximately  gaussian  with  mean  gq  and  variance 


s  )  =  &: 
o 


"'“i  •  "o> 


It  might  be  noted  that  this  variance  is  identical  with  the  one  that 


would  have  been  obtained  by  straightforward  application  of  naximura- 


likelihood  estimation  of  s  . 


The  predicted  error  at  time  t  =  t  Is  now  obtained  by  direct 


application  of  methods  used  in  Appendix  B.  Since  only  a  single  i^aromcter 


t=t 

c 


It  is  again  convenient  to  convert  the  summation  into  on  integral. 
Also,  it  is  convenient  to  let  t  =0.  Then 


and 


SR*"  = 
c 


dR(t-s)1 

OB  S=S 

■h=t 


a  T 


o  fh 

c  0 


bo) 


R'(t-s  )  dt 

o  J 


SPECIFIC  EXAMPLE 


As  in  Appendix  B,  ve  now  assume  that  the  form  of  R  is  given  by 


R(t-s)  =  +  ar^t-s) 


^  |  l-aru(t-B)  logj  l-or,((t-s)  +  ari((t-c 


(5D 


where  the  cr's  are  defined  as  in  Eq.  (26),  Appendix  D.  Then 


R'(t~s)  =  ~  =  -  a2  «-  log  l-ar4(t-s) 


(53) 


+  Li  (l+U  )  lOg^l  +  U  )  -  (l-u)log  (1-U) 

7Jt  L  o  0 


where 


U  =  O'),  s 

o  o 

UC  =  %^c  '  2o) 

U  -  6o) 

For  small  values  of  u  and  uq  it  lo  convenient,  to  replace  1ob(1+uq) 
and  log(l-u)  by  a  series.  A  particularly-sinplc  expression  results 

when  o',  =  0: 


C . 

ff  T 


W(1  - 


uc> 


u 

bO 


c,  r, 

(u^nr )+  .  • 
v  0 


(5M‘ 


From  this  it  appears  that  for  small 


5R"  decreases  as 


In  general, 


5R 

c  ■ 

G~  r 


r  k  -  log(l  -  u 


- -  3  ,  3  44 

,,  u  +u  u  -u 

2  .,22  o  ,  o  ^ 

v  li f  k(u  -u0,”r  ~~£~  ■■■) 


u*-u*  Ills' -ir) 


ijQ 


(55) 
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where  k  =  aja^ ■  For  positive  k  this  result  is  usually  smaller  than 
for  k  =  0.  Hence  curves  are  plotted  only  for  c*2  =  0.  These  are 
shown  in  Fig.  4. 
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Appendix  D 

PREDICTION  OF  SLBM  IN-PLANE  MOTION  DURING  BOOST  WHEN  PROPULSION 
to  vrsrtrH  mn  mNSTANT-THR'JST  AND  WHEN  ONLY  JTIME_  ORIGIN  AND 
~  INITIAL  POSITION  ARE  UNKNOWN 

The  results  of  Appendix  C  assume  complete  Knowledge  of  all  the 

parameters.  It  is  more  realistic  to  assume  that  the  initial  position 

parameter,  (y,  is  only  partially  Known  --  that  is,  one  docs  not  Know 

in  advance  when  the  SLBM  was  launched  nor  its  position  at  launch. 

In  this  section  it  is  therefore  assumed  that  a o/^,  and  a  (  are  Known 

exactly;  that  the  a  priori  probability  density  of  y  is  Russian  with 

a  mean  o  and  variance  «*;  and  that  s,  the  time  origin,  has  a  uniform 
m  ot 

distribution  as  in  Appendix  C.  [It  can  be  easily  demonstrated  that 
the  exact  form  of  the  a  priori  probability  distribution  of  y  or  of 
s  is  not  Important  after  the  first  few  seconds  of  observation  time.) 

Under  these  conditions,  the  ground  rar.i.e  o.  *ne  urn  •  missile 
can  be  expressed  in  the  fom: 

R(t)  =-  o x  '  f(t  -  s) 

where  f ( t  -  s)  is  an  exactly  Known  function.  The  a  posteriori  probu 
bility  of  y  and  s,  given  the  observed  values  oi  R,  is  then  g'-ven  »y 


p (s,  a  |R)  K  ;>(s,  ay 


t;(s.  y) 


where  k  is  again  a  normalising  constant,  and  where 


«(s,  oy)  - 


Ri  -  “i  -  fi'Li 


i  1 


V 


k6 

2 

The  noise  is  white,  gaussian,  with  variance  o',  as  in  Appendix  C,  and 
the  R^  are  the  observed  values  of  R(t).  If  the  noise  is  n(t)  and  if 
the  true  values  of  and  s  are  or  and  sq  respectively,  then  the 
observed  value  of  R  is 

R(t.)  =  Ri  =  or10  +  f(ti  -  sq)  +  n(tt) 


i 


and  therefore 


g(s,o1)  = 


1  V 


2  a 


2  )  aio  +  fl'Vs Q)  +  (50) 


i=l 


As  in  Appendix  C  we  expand  g(s,a  )  about  s  and  o  ,  and  vo  aait  terns 

O  i-J 

containing  n(t)  in  accordance  with  Woodward's  argument.  Under 
these  conditions  it  is  easily  shown  that 


g(Go'a10)  ° 


(59 


%  -  HSj  o 

00,  S  -S  cc  R--S 
1  c 


0 

°T  710  "l^lO 


(to) 


.2 

a  g- 

2. 

OO,  S  S 
1  0 

"l^lO 


n 

G. 

a 


since  n  -  t  /t 


(ol) 


.  2  ,  ’  J  4 

^ 6  '  _  n  t  /  a.  N 

2  ~  *  2  f  vti  ‘  G0) 

os  s-s  a  ,  , 

o  i-1 

*1  s10 


(te) 


1+7 


«. 


■  4  I  ''(‘i  -  ■.) 


(63) 


1=1  ■ 


“l^io 


-  s) 

where  f'v^  -  sq)  =  - s - 


8=8. 


Hence 


6(8,0.)  « - p 

1  2  a 


n 

r  <<yV2  - 2  ) 


1=1 


+  X  f'(trBo)  (B-Do)£ 


1  =  1 


(bl+) 


Also  In  view  of  the  assumptions  made  about  the  a  priori  probability 
distributions  of  or^  and  s,  we  have 


.  _L 

p(s,cr1)  =  ^  e  2aa 


(“iV2 


m 


By  means  of  seme  algebra  the  a  posteriori  distribution  p(d,q'1'R)  can 
then  be  put  In  the  form 


S^B,  Qfj) 

p(s,»1|R)  =  e 


(66) 


where 


61(6#Qf1) 


2.  2  \  2 

aio  Vi :  v  t 

*1*2.  2 

a  t.  +  o  t 

cr  1 


n  2 

-2  )  f,(t1-6o)(ol-Qf10)(B-Bo)  +yj[f'(t1'B0)]  (B-B0)  >(67) 

i=i  i=i  J 


The  mean  value  of  in  this  expreBBion  1b  a  linear  combination  of 
the  a  priori  mean  value  and  the  true  value  ot^q.  a  priori  mean 

value  would  be  different  from  the  true  value  only  If  the  measurement 
of  cr1  that  re suite  in  the  a  priori  distribution  were  biased.  There 
is  no  particular  reason  for  assuming  the  existence  of  Buch  a  bias 
here,  and  we  therefore  assume  that  or^  =  or^.  It  might  be  noted  that 

if  there  were  a  bias  in  the  measurement,  its  effect  would  in  any  case 

2  2  2 

become  negligible  if  a  t1  >>  a  r.  Unless  is  very  much  smaller 

than  a2 --in  which  case  one  would  be  justified  In  considering  ar^  to  be 
known  exactly--the  effect  of  the  bias  is  therefore  negligible  after 
only  a  few  seconds  of  observation  time. 

It  is  also  evident  that  if  aa  **  a 

2 

1  +  ~~  w  l  for  t  >  >  t  (66) 

a  1 


Hence,  if  t  is  on  the  order  of  .01  sec  or  less,  after  about  1  sec 
gjU,^)  is  given  to  a  very  good  approximation  by 


=  gU,^) 


i.e.,  the  a  posteriori  distribution  of  a  and  becomes  completely 
independent  of  a  priori  knowledge  of  Since  Eq.  (69)  also  applies 
to  the  case  of  no  a  priori  information  about  or^  which  represents  a 
least  favorable  situation,  this  expression  will  be  used  in  all  subse¬ 
quent  computations. 

The  variances  and  covariances  of  and  s  can  now  be  obtained 
direct^  from  Eq.  (69).  Note  that  this  equation  is  in  the  form: 


5l(B,c*i)  =  -  \  LA(<vV2  '  2B(Qf-Qrio)(s'So)  +  C(G'So)  j 


Hence 


("i  -  alOY 


2 


AC  -  B 


(B  - 


AC  -  B 


(70) 


(a. 


»10)<S 


.)  ■ 


B 


AC  -  B 


As  in  Appendix  C  the  prediction  error  at  time  t  is  obtained  from 

c 


r5R  1  - 


cJ 


|2 

cR  " 

a 

i _ 

(ori-Q-io)2  + 

c 

os 

_ 

12 


(s-s  )‘ 
0 


+  2 


dR 


00, 


oR 

c 

3i~ 


(orQ'io)(G'20) 


(7D 


where  the  partial  derivatives  are  all  evaluated  for  s  =  s  ,  a  -  , 

o’  10' 

t  =  t  .  If  R(t  -  s)  is  given  by  Eq.  ( 51) ,  Appendix  C.  then 


oR 


CiQ\ 


-  1 


dR 


os 


:S  =  -  1  10g(l-U  ) 

r\ 


(72) 


where  u  or,  ( t  -  s  ) . 
c  •»  c  o' 

The  two  summations  appearing  in  Eq.  (6 (j)  can  be  approximated  b> 


integrals,  as  usual.  Tiie  final  result  for  r5R  r'  Is  then- 

c 


a  a. 


“p1!  -  2  „  (or, ,)!>,( -u  )-g(u)]  +  -j.  [h(-u  )-hl  u) 1 
l  ofj  j  d  o  or.,  o 


[6R  ]2  +  2j  o2-o3  log(l-uc)  -o2ti-  ^  g(-uQ)-g(u)  +  Of2-of,  log(l-u  ) 


2 

y3 


"2  t(u+u0)[h(-u0)  '  h(u)]  -  2(u+uo)  g(-uQ)  -  g(u)  -  g(-u  )-g(u) 


2 

0  T 


This  expression  is  again  considerably  simplified  if  a  -  0,  and  this 
is  the  only  case  for  which  detailed  computations  have  been  performed 
In  this  case  [8R  is  independent  of  or  . 


The  results  for  this  case  calculated  by  a  digital  computer  from 


Eq.  (73)  for  several  values  of  s  are  presented  in  Fig.  3* 
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Appendix  E 
TYPICAL  TRAJECTORIES 

In  all  the  computations  performed  in  this  Memorandum,  a  typical 
SLBM  trajectory  has  been  assumed  when  specific  parameter  values  were 
required.  The  ground-range  function  corresponding  to  this  trajectory 
is  shown  in  Fig.  8.  The  equation  of  this  function  is 


R(t)  =  (1  -  .0143  t)  log(l  -  .0143  t)  +  .0143  tj  ( 70) 


Figure  8  also  presents  an  approximation  to  Eq.  (75)  by  means  of  a 
third-order  polynomial  function.  The  coefficients  of  this  polynomial 
have  been  chosen  so  that  the  polynomial  coincides  with  Eq.  (75)  for 
t  =  0,  20,  LO ,  and  00  sec.  The  equation  for  the  polynomial  approxi¬ 
mation  is: 


R(t)  =  120.66  l  +  4.75  tJ  +  .2558  tj 


ctj=  1970  ft/sec 


a*;. 0143  sec 


50,000 


R<  t )  =  120.66  1+4751  +  25581 
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Fig.  8  — Range  functions  for  typical  trajectory 
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